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Abstract 

After extensively reviewing general relativistic gravitomagnetism, 
both historically and phenomenologically, we review in detail the so- 
called magnetic components of gravitational waves (GWs), which have 
to be taken into account in the context of the total response functions 
of interferometers for GWs propagating from arbitrary directions. Fol- 
lowing the more recent approaches of this important issue, the analysis 
of such magnetic components will be reviewed in both of standard Gen- 
eral Theory of Relativity (GTR) and Scalar Tensor Gravity. Thus, we 
show in detail that such a magnetic component becomes particularly 
important in the high-frequency portion of the range of ground based 
interferometers for GWs which arises from the two different theories of 
gravity. Our reviewed results show that if one neglects the magnetic 
contribution to the gravitational field of a GW, approximately 15% of 
the potential observable signal could, in principle, be lost. 

Key words: Experimental studies of gravity Experimental tests of grav- 
itational theories Gravitational waves 
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1 General overview of gravitomagnetism 

The term "gravitomagnetism" [H [21 [3] (GM) commonly indicates the collec- 
tion of those gravitational phenomena regarding orbiting test particles, pre- 
cessing gyroscopes, moving clocks and atoms and propagating electromag- 
netic waves [H [5] which, in the framework of the Einstein's General Theory 
of Relativity [6] (GTR), arise from non-static distributions of matter and 
energy. In the weak-field and slow motion approximation, the Einstein field 
equations ^ of GTR, which is a highly non-linear Lorentz-covariant tensor 
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theory of gravitation, get linearized [8], thus looking like the Maxwellian 
equations of electromagntism [9]. As a consequence, a "gravitomagnetic" 
field Bg, induced by the off-diagonal components goi, i = 1,2,3 of the 
space-time metric tensor related to mass-energy currents, arises. In par- 
ticular, far from a localized rotating body with angular momentum S the 
gravitomagnetic field can be written as [TO] 

S-3(^S-r^f , (1) 

where G is the Newtonian gravitational constant and c is the speed of light 
in vacuum. It affects, e.g., a test particle moving with velocity v with a 
non-central acceleration [TO] 

= X Bg, (2) 

which is the cause of two of the most famous and empirically investigated 
GM effects, as we will see in Section 1.1 and Section 1.2. 




1.1 Historical overview 

The formal analogies between gravitation and electromagnetism date back 
to the early days of the Coulomb's force law |1H \T2\ [TO] between two non- 
moving pointlike electric charges (1785). Indeed, it follows faithfully the 
Newtonian inverse-square force law of gravitation [14j between two pointlike 
masses (1687) whose state of motion is, instead, irrelevant for its validity 
within the framework of classical mechanics. After the electrodynamical 
forces between current elements were discovered in 1820-1825 by Ampere 
[15] . the situation was reversed. Indeed, in 1870 Holzmiiller jtlGj, in order 
to study the motion of a test particle attracted by a fixed center, asked 
whether Newton's law might not be modified in much the same way as that 
in which Weber [n\ [TO] in 1846 had modified Coulomb's law for electric 
charges in an action-at-a-distance fashion by introducing velocity-dependent 
terms. He found that the trajectory is no longer closed, but can be described 
by a slowly precessing ellipse. Later, Tisserand [El |20] used this method, 
yielding a further, "magnetic" -like component of the Sun's gravitational field 
acting on the solar system's planets, to attempt-unsuccessfully-to explain 
the anomalous prograde perihelion precession of Mercur}^ of 43.98 arcsec 

^It found a natural explanation in 1915 by Einstein [21] within his GTR. Note that 
the static, "gravitoelectric" part of the Sun's field is required for the explanation of such 
a phenomenon: no mass-energy currents are involved. 
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cty~^, discovered in 1859 by Le Verrier [22]. Also Levy [23] worked in the 
same conceptual framework to solve that astronomical problem, but without 
success. The advent of Maxwell's electromagnetism [21], which is a linear 
field theory replacing the previous action-at-a-distance theorie^, did not 
discourage further attempts to use analogies of electromagnetic equations 
to solve gravitational problems. Maxwell himself |24j considered whether 
Newtonian gravity could be described by a vector field theory, but he did 
not succeed because of issues encountered with gravitational energy. A later 
attempt was due to Heaviside [26] in 1894 with his Maxwellian vector field 
theory of gravity. Among other things, he derived a planetary precession 
induced by the rotating Sun's GM dipole, but it was too small in magnitude 
and retrograde with respect to the Mercury's precession observed by Le 
Verrier. 

Another line of reasoning which yielded to consider gravitational forces 
induced by moving masses was that connected to the need of explaining 
the origin of inertia. In particular, the idea that rotating bodies may ex- 
ert not only the static Newtonian gravitational force but also an additional 
"frame-dragging" on test particles was probably due to Mach [27]. He spec- 
ulated that even for relative rotations centrifugal forces arise due to some, 
unspecified gravitational interaction with the masses of the Earth and of 
the other celestial bodies [28], certainly quite larger than the mass of the 
famous Newton's bucket in relative rotation with respect to the water in- 
side. Later, in 1896 the Friedlander brothers [29] expressed the conviction 
that the properties of inertia and gravitation should be finally derived from 
a unified law. A Maxwell-like theory of gravitation, proposed to explain the 
origin of inertia in the framework of the Mach's principle, is due to Sciama 

m- 

After the birth of the Einstein's Special Theory of Relativity (STR) in 
1905 [31], the problem of a "magnetic" -type component of the gravitational 
field of non-static mass distributions was tackled in the framework of the 
search for a consistent relativistic theory of gravitation [32]. Indeed, bringing 
together Newtonian gravitation and Lorentz invariance in a consistent field- 
theoretic framework necessarily requires the introduction of a "magnetic" - 
type gravitational field of some form [331 EH [35] . 

With a preliminary and still incorrect version of GTR, Einstein and Besso 
in 1913 [36j calculated the node precession of planets in the field of the ro- 
tating Sun; the figures they obtained for Mercury and Venus were incorrect 

•^A recent action-at-a-distance gravity tiieory, based on the scalar velocity-dependent 
Weber-type potential, is due to Assis [25| . 
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also because they used a wrong value for the solar mass. Soon after GTR 
was put forth by Einstein (1915) [37], de Sitter [38] in 1916 used it to prelim- 
inarily work out the effects of Sun's rotation on planets' perihelia, although 
he restricted himself to ecliptic orbits only; his result for Mercury (—0.01 
arcsec cty~^) was too large by one order of magnitude because he assumed a 
homogenous and uniformly rotating Sun. In 1918 Thirring [39] analyzed in 
a short article the formal analogies between the Maxwell equations and the 
linearized Einstein equations. Later [391 140j . Thirring computed the cen- 
trifugal and Coriolis-like gravitomagnetic forces occurring inside a rotating 
massive shell. Lense and Thirring [H] in 1918 worked out the gravitomag- 
netic effects on the orbital motions of test particles outside a slowly rotating 
mass; in particular, they computed the gravitomagnetic rates of both the 
satellites of Mars (Phobos and Deimos), and of some of the moons of the 
giant gaseous planets. They found for the longitude of the ascending node 
O a prograde precession, while for the argument of pericenter uj a retrograde 
precession occurs 



where a, e, I are the semimajor axis, the eccentricity and the inclination of 
the test particle's orbital plane to the central body's equator, respectively. 

Another well known general relativistic gravitomagnetic effect consists 
of the precession of a gyroscope moving in the field of a slowly rotating body. 
It was derived in 1959 by Pugh [35] and in 1960 by Schiff [Ml iZl HH] . The 
possible Machian character of the Schiff effect was discussed by Rindler [49], 
and Bondi and Samuel [50j. 

Cosmological GM and Mach's principle have been recently treated by 
Schmid [SB [22]. 

Certain subtle issues concerning the gravitomagnetic effects inside a ro- 
tating massive shell were solved by Pfister and Braun in 1985 |S3]. 

For seeming analogies between Maxwellian electromagnetism and the 
fully non-linear equations of GTR, see the works by Matte |54j and, more 
recently, by Costa and Herdeiro [55]- In this framework, Pascual-Sanchez 
[56^1 discussed the non-existence of a GM dynamo in the linearized, weak-field 

^However, in August 1917 Einstein j42J wrote to Thirring that he calculated the 
Coriolis-type field of the rotating Earth and Sun, and its influence on the orbital ele- 
ments of planets (and moons). A detailed history of the formulation of the so-called 
Lense- Thirring effect has recently been outlined by Pfister [43] ; according to him, it would 
be more fair to speak about an Einstein-Thirring-Lense effect. 
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and slow-motion approximation of GTR. Tartaglia and Ruggiero [57] inves- 
tigated the possible occurrence of a GM analog of the Meissner effect using 
the same approximation for GTR. The non-existence of such a phenomenon 
in gravitation and of other putative GM effects has been demonstrated by 
Pascual- Sanchez [55] . 

1.2 Experimental/observational overview 

About empirical investigations of possible gravitational effects due to moving 
bodies, in 1896 I. Friedlander [29] was the first to perform an experiment 
concerning a putative gravitational influence of moving bodies. He used as a 
source a rapidly rotating heavy fly wheel and tried-unsuccessfully-to detect 
its gravitational effects on a torsion balance mounted above the fly wheel, in 
line with its axis. Later, in 1904 Foppl |59] looked for possible Coriolis-like 
gravitational dragging effects induced on a gyroscope made of two heavy fly 
wheels by the whole rotating Earth as a source. 

Moving to more recent epochs, soon after the dawn of the space age with 
the launch of Sputnik in 1957 it was proposed by Soviet scientists to directly 
test the general relativistic Lense-Thirring effect with artificial satellites or- 
biting the Earth. In particular, V.L. Ginzburg [60l [HI [62] proposed to use 
the perigee of a terrestrial spacecraft in highly elliptic orbit, while A.F. Bo- 
gorodskii j63j considered also the node. In 1959 Yilmaz |64j . aware of the 
aliasing effect of the much larger classical precessions induced by the non- 
sphericity of the Earth, proposed to launch a satellite in a polar orbit to 
cancel them. About twenty years later, in 1976 van Patten and Everitt 
[65| |66] suggested to use a pair of drag-free, counter-orbiting terrestrial 
spacecraft in nearly polar orbits to detect their combined Lense-Thirring 
node precessions. In 1977-1978 Cugusi and Proverbio [671 EH] suggested to 
use the passive geodetic satellite LAGEOS, in orbit around the Earth since 
1976 and tracked with the Satellite Laser Ranging (SLR) technique, along 
with the other existing laser-ranged targets to measure the Lense-Thirring 
node precession. In 1986 Ciufolini |69j proposed a somewhat simpler ver- 
sion of the van Patten-Everitt mission consisting of looking at the sum of the 
nodes of LAGEOS and of another SLR satellite to be launched in the same 
orbit, apart from the inclination which should be switched by 180 deg in 
order to minimize the competing classical precessions due to the centrifugal 
oblateness of the Earth. lorio |70] showed that such an orbital configuration 
would allow, in principle, to use the difference of the perigees as well. Tests 
have started to be effectively performed later by Ciufolini and coworkers 
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[71] with the LAGEOS and LAGEOS II satehitefl according to a strat- 
egy by CiufoUni |72) involving the use of a suitable linear combination of 
the nodes of both satellites and the perigee uj of LAGEOS II in order to 
remove the impact of the first two multipoles of the non-spherical gravita- 
tional potential of the Earth. Latest tests have been reported by Ciufolini 
and Pavlis [73^ , Lucchesi [75] and Ries and coworkers [76] with only the 
nodes of both the satellites according to a combination of them explicitly 
proposed by lorio [77j . The total uncertainty reached is still matter of de- 
bate [TSl ESI EQl ED Ell [831 El] because of the lingering uncertainties in the 
Earth's multipoles and in how to evaluate their biasing impact; it may be as 
large as ^ 20 — 30% according to conservative evaluations [78 | [8T | [82 | [83 | [84] . 
while more optimistic views [73\ I7H [76] point towards 10 — 15%. Several 
authors [HI ESJ [571 ESl [89] explored the possibility of using other currently 
orbiting SLR geodetic satellites in addition to LAGEOS and LAGEOS II. 
A new SLR geodetic satellite, named LARES, should be launched by the 
Italian Space Agency (ASI) in 2010 to improve the present-day accuracy of 
the Lense-Thirring tests by combining its node with those of the existing 
LAGEOS and LAGEOS II |9l]. The claimed accuracy is 1% [90], but also in 
this case the realistic level of uncertainty may be quite larger |92[ [93| [Ml [95] 
because of the relatively low orbit of LARES with respect to LAGEOS and 
LAGEOS II which should bring into play the systematic alias by several 
non-perfectly known Earth's multipoles. 

In 2006 a preliminary test in the gravitational field of Mars with the 
Mars Global Surveyor (MGS) has been performed by lorio [96l [97]. He 
interpreted certain features of the time series of the out-of-plane portion N 
of the MGS orbit involving its node in terms of the Lense-Thirring effect. 
In particular, the average of the Root-Mean-Square (RMS) orbit-overlap 
differences of the out-of-plane component of the MGS path over 5 years 
agrees with the predicted average Lense-Thirring out-of-plane shift over the 
same time span within a few percent, while a linear fit to the complete 

time series for the entire MGS data set shows an agreement with the 
corresponding predicted Lense-Thirring signal at a ~ 40%. A debate about 
the validity of such an interpretation arose [98J , and it is still ongoing [99] . 
The possibility of designing a dedicated mission to Mars has been recently 
considered by lorio [100] . 

lorio and Lainey [101] revisited the original proposal by Lense and Thirring 
[44] concerning the system of Jovian moons in view of recent advances in 

^LAGEOS II was launched in 1992, but its orbital configuration is different from that 
proposed in [69 . . 
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orbit determination of the four large Galilean satellites. The possibility of 
using fruitfully them seem still to be premature. 

Concerning the Sun's GM field and the inner planets, the situation is 
nowadays more favorable than in the past [381 168j . Indeed, the astronomer 
Pitjeva |lU2tllU3] has recently fitted the full set of dynamical force models of 
the planetary motions of the EPM ephemerides to about one century of data 
of several types by estimating, among other things, corrections Aro to the 
standard Newtonian/Einsteinian secular precessions of the longitude of the 
perihelia of all the rocky planets. In doing so she did not model the solar GM 
field, so that such corrections, by construction, are well suited to test the 
Lense-Thirring effect [1041 11051 1106j . The magnitude of the predicted Lense- 
Thirring perihelion precessions, although one order of magnitude smaller 
than what argued in earlier studies \38\ I68j . is about of the same order 
of magnitude of, or even larger than, the present-day uncertainty in the 
estimated Aw. In particular, it has been noted |106j that the Lense-Thirring 
perihelion precession of Venus amounts to vjut = —0.0003 arcsec cty~^, 
while the estimated correction for Venus is Azu = —0.0004 it 0.0001, where 
the quoted uncertainty is the 1 — a statistical error. Thus, the existence 
of the Lense-Thirring effect would be confirmed at a 25% level, although 
caution is in order because the realistic uncertainty in Adj might be up 
to 5 times larger. The systematic bias due to the mismodelling of other 
competing classical effects would be less relevant than in the LAGEOS- 
LAGEOS II case. The proposed space-based Astrodynamical Space Test 
of Relativity using Optical Devices (ASTROD) mission [107j has, among 
its scientific goals, also the accurate determination of the Sun's angular 
momentum through a ~ 1% measurement of the gravitomagnetic time delay 
on electromagnetic waves. 

Soon after the formulation of the Schiff effect, in 1961 Fairbank and Schiff 
|108j submitted to NASA a proposal for a dedicated space-based experiment 
aimed to directly measure it. Such an extremely complicated mission, later 
named Gravity Probe B (GP-B) |109t 1110] . consisted of a drag-free, liquid 
helium-cooled spacecraft moving in a polar, low orbit around the Earth and 
carrying onboard four superconducting gyroscopes whose GM precessions 
should have been detected by Superconducting Quantum Interference De- 
vices (SQUID) with an expected accuracy of 1% or better. It took 43 years 
to be implemented since GP-B was finally launched on 20 April 2004; the 
science data collection lasted from 27 August 2004 to 29 September 2005, 
while the data analysis is still ongoing |11H I112j . It seems that the final 
accuracy obtainable will be not so good as initially hoped because of the 
occurrence of unexpected systematic errors [1131 111^ I115j . At present, GP- 
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B team report^ a statistical error of ~ 14% and systematic uncertainty of 
~ 10%. In 1975 Haas and Ross |116j proposed to measure the angular mo- 
menta of the Sun and Jupiter by exploiting the Schiff effect with dedicated 
spacecraft-based missions, but such a proposal was not carried out so far. 

All the previously reviewed attempts aim to obtain direct tests of some 
GM effects. However, according to Nordtvedt [1171 1118) . the GM interac- 
tion would have already been observed, with a relative accuracy of 1 part 
to 1000, in comprehensive fits of the motions of several astronomical and 
astrophysical bodies like satellites, binary pulsars and the Moon. In fact, 
Nordtvedt does not refer to the effects considered so far, caused by the ro- 
tation of the body which acts as source of the gravitational field ( "intrinsic" 
GM). Instead, he primarily deals with some GM long-periodic harmonic 
perturbations affecting the Earth-Moon range induced by the translational 
GM mass currents due to the orbital motion of the Earth-Moon system 
around the Sun ("extrinsic" GM); Lunar Laser Ranging (LLR) would have 
measured them with a 0.1% accuracy |119| . However, Kopeikin argues that 
LLR would not be able to detect genuine GM signatures which are not spu- 
rious, gauge-dependent effects |120] . For other works about such an issue, 
see \121\ \122\ I123j . Concerning the possibility of directly measuring the 
Lense-Thirring precessions of the Moon's motion due to the Earth's angular 
momentum, it has been recently proven to be still unfeasible by lorio |124| 
because of the too large level of uncertainty in several competing classical 
effects. A test of extrinsic gravitomagnetism concerning the deflection of 
electromagnetic waves by Jupiter in its orbital motion has been reported in 
a dedicated analysis of radiointerferometric data by Fomalont and Kopeikin 
|125j . but also such a test is controversial: see the WEB page by Will at 
|http:/ /physics. wustl.edu/cmw/SpeedofGravity.htmT, 

For other proposals to directly detect various aspects of the (intrinsic) 
GM field in Earth-based laboratory and space-based experiments, see, e.g., 
the book by lorio |126j . Extensive overviews of the importance of GM 
in astrophysical contexts like accretion disks around compact objects and 
relativistic jets in quasars and galactic nuclei can be found, e.g., in the book 
by Thorne, Price and MacDonald jl27!, and in Section E of the book by 
Ruffini and Sigismondi |128j . 



See on the WEB: http://einstein.stanford.edu/| 
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2 Magnetic component in the gravitational field 
of a gravitational wave 

2.1 The importance of gravitational waves: a new window 
into the Universe 

The data analysis of interferometric Gravitational Waves (GWs) detectors 
has nowadays been started, and the scientific community hopes in a first 
direct detection of GWs in next years; for the current status of GWs in- 
terferometers see Ref. [130j . In such a way, the indirect evidence of the 
existence of GWs by Hulse and Taylor [131] . Nobel Prize winners, will be 
confirmed. Detectors for GWs will be important for a better knowledge of 
the Universe [130j and also because the interferometric GWs detection will 
be the definitive test for GTR or, alternatively, a strong endorsement for 
Extended Theories of Gravity |132j . In fact, if advanced projects on the 
detection of GWs improve their sensitivity, allowing the Scientific Commu- 
nity to perform a GW astronomy, accurate angle- and frequency-dependent 
response functions of interferometers for GWs arising from various theories 
of gravity will permit to discriminate among GTR and extended theories of 
gravity. This ultimate test will work because standard GTR admits only 
two polarizations for GWs, while in all extended theories the polarizations 
are, at least, three, see |132j for details. 

On the other hand, the discovery of GW emission by the compact binary 
system composed by two Neutron Stars PSR1913+16 |131j has been, for 
physicists working in this field, the ultimate thrust allowing to reach the 
extremely sophisticated technology needed for investigating in this field of 
research [130| . 

Gravitational Waves are a consequence of Einstein's GTR |133j . which 
presuppose GWs to be ripples in the space-time curvature travelling at light 
speed [13411135] . Only asymmetric astrophysics sources can emit GWs. The 
most efficient are coalescing binaries systems, while a single rotating pulsar 
can rely only on spherical asymmetries, usually very small. Supernovae 
could have relevant asymmetries, being potential sources |130j . 

The most important cosmological source of GWs is, in principle, the 
so called stochastic background of GWs which, together with the Cosmic 
Background Radiation (CBR), would carry, if detected, a huge amount of 
information on the early stages of the Universe evolution \136\ 11371 I138j . 
The existence of a relic stochastic background of GWs is a consequence of 
generals assumptions. Essentially it derives from a mixing between basic 
principles of classical theories of gravity and of quantum field theory. The 
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strong variations of the gravitational field in the early universe amplify the 
zero-point quantum oscillations and produce relic GWs. It is well known 
that the detection of relic GWs is the only way to learn about the evolution 
of the very early universe, up to the bounds of the Planck epoch and the 
initial singularity [1391 114U| . It is very important to stress the unavoidable 
and fundamental character of this mechanism. The model derives from the 
inflationary scenario for the early universe [141] . which is tuned in a good 
way with the WMAP data on the CBR (in particular exponential inflation 
and spectral index ~ 1 [142j ). Inflationary models are cosmological models 
in which the Universe undergoes a brief phase of a very rapid expansion in 
early times |141j . In this context the expansion could be power-law or expo- 
nential in time. Such models provide solutions to the horizon and flatness 
problems and contain a mechanism which creates perturbations in all fields 
[1391 11401 1141j . Important for our case is that this mechanism also provides 
a distinctive spectrum of relic GWs [1391 1140j . The GWs perturbations arise 
from the uncertainty principle and the spectrum of relic GWs is generated 
from the adiabatically-amplified zero-point fluctuations \139\ 1140] . 

Regarding the potential GW detection, let us recall some historical notes. 

In 1957, F.A.E. Pirani, who was a member of the Bondi's research group, 
proposed the geodesic deviation equation as a tool for designing a practical 
GW detector [Tii] . 

In 1959, Joseph Weber studied a detector that, in principle, might be 
able to measure displacements smaller than the size of the nucleus |145j . He 
developed an experiment using a large suspended bar of aluminum, with a 
high resonant Q at a frequency of about 1 kHz. Then, in 1960, he tried 
to test the general relativistic prediction of gravitational waves from strong 
gravity collisions [146] and, in 1969, he claimed evidence for observation of 
gravitational waves (based on coincident signals) from two bars separated by 
1000 km [147j . He also proposed the idea of doing an experiment to detect 
gravitational waves using laser interferometers [147j . In fact, all the modern 
detectors can be considered like being originated from early Weber's ideas 

m- 

At the present time, in the world there are five cryogenic bar detec- 
tors have been built to work at very low temperatures (< 4K): Explorer 
at CERN, Nautilus at Frascati INFN National Laboratory, Auriga at Leg- 
naro National Laboratory, Allegro at Luisiana State University and Niobe 
in Perth [130\ . Instrumental details can be found in [130j and references 
within. Spherical detectors are the Mario Schenberg, which has been built 
in San Paolo (Brazil) and the MiniGRAIL, which has been built at the 
Kamerlingh Onnes Laboratory of Leiden University, see |130j and references 
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within. Spherical detectors are important for the potential detection of the 
scalar component of GWs that is admitted by Extended Theories of Gravity 
[148) . In the case of interferometric detectors, free falling masses are inter- 
ferometer mirrors which can be separated by kilometres (3km for Virgo, 
4km for LIGO). In this way, GW tidal force is, in principle, several order 
of magnitude larger than in bar detectors. Interferometers have very large 
bandwidth (10-10000 Hz) because mirrors are suspended to pendulums hav- 
ing resonance in the Hz region. Thus, above such a resonance frequency, 
mirrors works, in a good approximation, like freely falling masses in the 
horizontal plane |130j . 

Recently, starting from the analysis in Ref. |149j . some papers in lit- 
erature have shown the importance of the gravitomagnetic effects in the 
framework of the GWs detection too [150^ 1151] I152j . In fact, the so-called 
magnetic components of GWs have to be taken into account in the context 
of the total response functions of interferometers for GWs propagating from 
arbitrary directions, [1491 11501 11511 [T52j . In next analysis we will show that 
such a magnetic component becomes particularly important in the high- 
frequency portion of the range of ground based interferometers for GWs 
which arises from standard GTR. 

In a recent paper, the magnetic component has been extended to GWs 
arising from scalar-tensor gravity too |153| . In particular, in Ref. |153| it 
has been shown that if one neglects the magnetic contribution considering 
only the low-frequency approximation of the electric contribution, a portion 
of about the 15% of the signal could be, in principle, lost in the case of 
Scalar Tensor Gravity too, in total analogy with the standard case of GTR 
[T49l[T50l[T5n[T52] . 

For the sake of completeness, such a case will be included in the following 
discussion. 

2.2 The magnetic component of GWs in standard GTR 

In a laboratory environment on Earth coordinate systems in which the space- 
time is locally flat are typically used, and the distance between any two 
points is given simply by the difference in their coordinates in the sense 
of Newtonian physics [1531 1154^ I155j . In this frame, called the frame of 
the local observer, GWs manifest them-self by exerting tidal forces on the 
masses (the mirror and the beam-splitter in the case of an interferometer 
[MKISIKIMIIIM], see Figure 1. 

The presence and importance of the so-called magnetic components of 
GWs in the framework of GTR was emphasized by Baskaran and Grishchuk 
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Figure 1: photons can be launched from the beam-sphtter to be bounced 
back by the mirror 



that computed the correspondent detector patterns |149j . while more de- 
tailed angular and frequency dependences of the response functions for the 
magnetic components have been given, with a specific application to the 
parameters of the LIGO and Virgo interferometers, in Refs. [150^ I15H 1152] , 
Thus, following Refs. [150^ I15H I152j . we will, now, work with G = 1, c = 1 
and h= 1 and can call h+{tu + zu) and {tu + zu) the weak perturbations 
due to the + and the x polarizations which are expressed in terms of syn- 
chronous coordinates tu, xtt,yu, zu in the transverse-traceless (TT) gauge. 
In this way, the most general GW propagating in the zu direction can be 
written in terms of a plane monochromatic wave [1501 11511 1152j 

h^u{ttt + ztt) = h+{ttt + ztt)eiiu' + /ix {ttt + ztt)ei^J = 

(4) 

= h+o exjpiuj{ttt + ztt)e'i^^ + /ixo exp iuj{ttt + ztt)e%\ 

and the corresponding line element will be 

ds^ = dtff- — dzu — (1 + hj^)dx1^ — (1 — hj^)dyj^ — 2hxdxttdxtt- (5) 

The wordlines xttiUtt-, zu = const, are timelike geodesies representing the 
histories of free test masses [1491 1150^ 11511 1152j . The coordinate transforma- 
tion x" = from the TT coordinates to the frame of the local observer 
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t = tu + \{xlt - ytt)h+ - \xttytthx 
x = xtt + \xtth+ - ^ytthx + \xttztth+ - \yttZtthy, 
y = ytt + \ytth+ - ^xuh^ + \yttztth+ - \xttztthx 

z = ztt- jixlt - ytt)h+ + ^xttytthx, 

where it is /i+ = and hx = The coefficients of this transformation 
(components of the metric and its first time derivative) are taken along 
the central wordline of the local observer [1491 11501 HSlj I152j . It is well 
known from Refs. |149tll5(Jlll51t[T^ that the linear and quadratic terms, as 
powers of xf^, are unambiguously determined by the conditions of the frame 
of the local observer, while the cubic and higher-order corrections are not 
determined by these conditions. Thus, at high-frequencies, the expansion in 
terms of higher-order corrections breaks down [149^ 1150^ 1151^ I152j . 

Considering a free mass riding on a timelike geodesic {x = h, y = I2, 
z = h) [1491 [IMI [ini [US], eq. dH) defines the motion of this mass with 
respect to the introduced frame of the local observer. In concrete terms one 
gets 

x{t) = li + \[lih+{t) - hhx {t)] + Ihkh+it) + Ihhhx (t) 

y{t) = h- khh+{t) + hhx it)] - ^hkh+it) + Ihhhx (t) (7) 

zit) = h-i[{ll- ll)h+{t) + 2hhhx (t), 

which are exactly eqs. (13) of Ref. [149| rewritten using the notation of 
Refs. |1501 I15H I152j . In absence of GWs, the position of the mass is 
(hyh^h)- The effect of the GW is to drive the mass to have oscillations. 
Thus, in general, from eq. ([7]) all three components of motion are present 
[T501 [T5H [T52] . Neglecting the terms with h+ and Ax in eq. the 
traditional equations for the mass motion are obtained [1491 11501 115H 1152] 

x{t) = h + ^[hh+{t) - hhxit)] 
y{t) = I2 - ^[hh+it) + hhxit)] (8) 
z{t) = k. 
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Clearly, this is the analogous of the electric component of motion in electro- 
dynamics [1491 11501 11511 1152j , while equations 



x{t) = li + l 



\lihh+{t) + IhhhAt) 



y{t) = h-h 



(9) 



z{t) = k- Mil - ll)h+{t) + 2hhh^{t) 



are the analogous of the magnetic component of motion. One could think 
that the presence of these magnetic components is a frame artefact due 
to the transformation eq. ([6]), but in Section 4 of Ref. |149j eq. ([7]) have 
been directly obtained from the geodesic deviation equation too, thus the 
magnetic components have a real physical significance. The fundamental 
point of Ref. [149t \\b[)\ \\b\\ I152j is that the magnetic components become 
important when the frequency of the wave increases but only in the low- 
frequency regime. This can be understood directly from eq. ([7]). In fact, 
using eq. @ and eq. ([6]), eq. ^ become 



x{t) = h + \[lih+{t) - hh-^it)] + \lihujh+{t - f ) + Ihhu^K {t - f ) 
y{t) = h- \[hh+{t) + /l/lx (t)] - \hhioh+{t - f ) + \hhujh-^{t - f ) 
z{t) = h-Wi- ll)ujh+{t - f ) + 2hl20jh^ it - f ). 



wavelength goes to infinity [149^ I150t I15H [T52] . while, at high-frequencies, 
the expansion in terms of ulilj corrections, with i,j = 1,2,3, breaks down 



Now, let us compute the total response functions of interferometers for 
the magnetic components in standard GTR. 

Equations eq. ([7]), that represent the coordinates of the mirror of the 
interferometer in presence of a GW in the frame of the local observer, can 
be rewritten for the pure magnetic component of the + polarization as 




[IMllinollEIllISS]. 



x{t) = li + \hhh+{t) 



y{t) = h-\ 



\hhh+{t) 



(11) 



z{t) = h-\{il-il)h+{t), 
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Figure 2: a GW propagating from an arbitrary direction 

where /i,/2 o-f^d h ^-re the unperturbed coordinates of the mirror. 

To compute the response functions for an arbitrary propagating direc- 
tion of the GW, one recalls that the arms of the interferometer are in general 
in the it and 17" directions, while the x,y, z frame is adapted to the prop- 
agating GW (i.e. the observer is assumed located in the position of the 
beam splitter). Then, a spatial rotation of the coordinate system has to be 
performed: 

u = —X cos 9 cos (t> + y sin (f) + z sin 9 cos cj) 

V = —X cos9 sincj) — y cos (j) + zsm9 sincf) (12) 

w = xsm9 + zcos9, 

or, in terms of the x,y,z frame: 

x = —u cos 9 cos (j) — V cos 9 sin (p + w sin 9 

y = n sin (/) — u cos (/> (13) 

z = usm9 coscj) + V sm9 sincj) + w cos9. 

In this way, the GW is propagating from an arbitrary direction 1^ to the 
interferometer (see Figure 2). As the mirror of eq. (|lip is situated in the 
u direction, using eq. (fTTI) . eq. (fT2]) and eq. (fT3]l the u coordinate of the 
mirror is given by 
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where 



u = L + ^L^Ah+{t), (14) 



A = sin 6 cos 0(cos^ 6 cos^ (j) — sin^ 0) (15) 



and L = y^/f + + ^3 is the length of the interferometer arms. 
The computation for the v arm is similar to the one above. Using eq. 
(jlip , eq. ()12p and eq. ()13p , the coordinate of the mirror in the v arm is 

V = L + ^L'^Bh+{t), (16) 

where 

B = sin 9 sin (cos^ cos^ (f> — sin^ (/>) . (17) 

Eq. (I14p and eq. (|16p represent the distance of the two mirrors of the 
interferometer from the beam-splitter in presence of the GW (note that only 
the contribution of the magnetic component of the + polarization of the GW 
is taken into account). They represent particular cases of the more general 
form given in eq. (33) of jl49j . 

A signal can also be defined in the time domain (T = L in our notation) 

Ml^^^}^,(A-B)hAt). (18) 

The quantity eq. (fT8|) can be computed in the frequency domain by 
using the Fourier transform of defined by 



h+{uj)= / dth+{t)exp{iujt), (19) 
J —00 

obtaining 



where the function 

H+{^) = -UuL{A-B) 



UojL sin e [(cos2 e + sin 2(/.i±^f^ )] (cos - sin , 



(20) 



is the total response function of the interferometer for the magnetic com- 
ponent of the + polarization [Ti9| [TSOl [T5T|[T52]. 
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Figure 3: the angular dependence of the response function of the LIGO 
interferometer to the magnetic component of the + polarization for / = 8000 
Hz 

In the above computation the theorem on the derivative of the Fourier 
transform has been used. 

The angular dependence of the response function eq. (j20p of the LIGO 
interferometer to the magnetic component of the + polarization for / = 8000 
Hz is shown in Figure 3. 

The analysis can be generalized for the magnetic component of the x 
polarization too. In this case, eq. ([7]) can be rewritten for the pure magnetic 
component of the x polarization as [1491 11501 115H I152j 

x{t) = li + ^hhhxit) 

y{t) = l2 + \hhhAt) (21) 

z{t) = h-\hl2h,,{t). 
Using eq. ([2T|) , eq. ([12]) and eq. , the u coordinate of the mirror in the 
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u arm of the interferometer is given by 



u = L + ^L^Ch^{t), (22) 

where 

C = —2 cos 6 cos^ (j) sin 6 sin (23) 

while the v coordinate of the mirror in the v arm of the interferometer is 
given by 

v = L+^L^DK{t), (24) 

with 

= 2 cos cos (/) sin sin^ </). (25) 

Thus, with an analysis similar to the one of previous Sections, it is possible 
to show that the response function of the interferometer for the magnetic 
component of the x polarization is [1491 11501 115H I152j 

^r^agn(^) = -i^T{C - D) = 

(26) 

= —iijoL sin 2(/)(cos (f) + sin </>) cos 9, 



The angular dependence of the response function eq. (j26l) of the LIGO 
interferometer to the magnetic component of the x polarization for / = 8000 
Hz is shown in Figure 4. 

From Figure 3 and Figure 4, it looks clear that if one neglects the mag- 
netic contribution, approximately 15% of currently observable signal could, 
in principle, be lost. 



2.3 The magnetic component of GWs in Scalar Tensor Grav- 

ity 

In the framework of Scalar Tensor Gravity, the TT gauge can be extended 
to a third polarization [Ml [ISSl UMl (ISZl EMI [EH HBO]. In this way, the 
total perturbation of a gravitational wave propagating in the z— direction 
in this gauge is [IMl [1531 1160J 

V(i + z) = A+{t + z)eW +A^{t + z)e^^J + ^{t + z)eW. (27) 

The term A^(i + z)e]ii) + A^{t + z)e^^^ describes the two standard 
(i.e. tensor) polarizations of gravitational waves which arises from General 
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Relativity in the TT gauge, see previous Subsection, while the term $(t + 
z)el^} is the extension of the TT gauge to the scalar-tensor case. 

For a purely scalar GW the metric perturbation eq. ([27]) reduces to 
[1321 1531 [160] 

V = Ht^ (28) 
and the correspondent line element is |132l I153[ I160j 

ds"^ = dt^ - dz^ - (1 + - (1 + (29) 

with $ = $oe*'^(*+'^^ 

Again, the wordlines x,y,z = const, are timelike geodesies representing 
the histories of free test masses, see the analogy with tensor waves in previous 
Subsection. In this case, the coordinate transformation = from 
the TT coordinates to the frame of the local observer is |153j 

t = ttt + \{xl - yl)^ 

X = Xtt + ^Xtt^ + ^XttZtt^ 

(30) 

y = ytt + \ytt^ + \yttztt^ 

z = ztt- \{xl-yl)^, 

where it is 6 = ^ , see previous Subsection and |153j . 

Now, if one considers a free mass riding on a timelike geodesic [x = li, 
y = I2, z = I3), eq. ([30[) define the motion of this mass due to the scalar 
GW with respect to the introduced frame of the local observer. Thus, one 
gets 

xit) = h + lh<^{t) + Ihh^it) 

y{t) = l2 + \hHt) + \hhHt) (31) 

z{t) = h-i^{ii-ii)m- 

In absence of scalar GWs the position of the mass is {h,h,h)- Again, 
the effect of the scalar GW is to drive the mass to have oscillations. Thus, 
in general, from eq. ([3ip all three components of motion are present. 

Neglecting the terms with $ in eq. ([3ip . the traditional equations for 
the mass motion due to the scalar GW are obtained [I53j 
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x{t) = /i + i/i$(i) 

y{t) =l2 + Ih-^it) (32) 
Z{t) = k. 

This is the analogous of the electric component of motion in electrodynamics 
(see previous Subsection), while equations 

x{t) =h + ^hk^t) 

y(t) = l2 + ^hh^it) (33) 

z{t) = h-i{ii-q)m, 

are the analogue of the magnetic component of motion. 

Thus, the magnetic component becomes important when the frequency 
of the wave increases in this case too, but only in the low-frequency regime 
in analogy with the standard tensor case. 

Even in this scalar case, one could think that the presence of this mag- 
netic component is a frame artefact due to the transformation eq. (|30p . 
but now we show that eq. (p3|) can be directly obtained from the geodesic 
deviation equation too, proving that the magnetic components have a real 
physical significance. 

Following ^153], let us focus the attention on the geodesic deviation ex- 
tended to second order approximation. The derivation of the geodesic de- 
viation equations is usually based on a two parameter family of timelike 
geodesic x"(r, r). In general and in the lowest approximation, the geodesic 
deviation equations are given by [1531 1160] 

^ = Rj^u^u-^nP. (34) 

The vector is the unit tangent vector to the geodesic and is the 
separation vector between two nearby geodesies 

u"{T,r) = |r=const., (35) 

n"'{T,r) = — — |t=const. • (36) 
or 

It is also assumed that the central geodesic line corresponds to r = 
while the second nearby geodesic corresponds to r = rg ^153j . Then, -R^,^^ 
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in eq. (|34p is the curvature tensor calculated along the central geodesic and 
^ the covariant derivative calculated along that line [153] . 

To discuss the magnetic component of motion in the field of a scalar 
GW, we need the geodesic deviation equations extended to the next ap- 
proximation. These equations have been obtained in |153] . Let us introduce 
the closely related vector w°' 

-" = ^ = ->'^ = ^ + r^y-^. (37) 

This vector obeys the equations |153j 



^ = Rj^u'^u^w^ + {Rj^.^-Rj.^^)u'^uPu^u^+ARjy^ (38) 
Defining the vector 

A^" = ron" + \rlw^, (39) 
eq. ([M|) and eq. ([37|) can be combined obtaining |153j 



= R^p^^u''u^NP+{Rj^.,-R^,l^p)u''u^N^N'+2Rjy^—N^+0{7 

(40) 

Thus, it is possible writing the expansion of x"(r, ro) in terms of N'^ 
[153] 

x°(r, ro) = a;"(r, 0) + iV° - T^N^N^ + 0(rl). (41) 

This formula shows that in the frame of the local observer (in which it 
is = along the central geodesic line |153j ) the spatial components of 
A'^" will directly give the time-dependent position of the nearby test mass. 
According to eq. (flO]l . these positions include the next-order corrections, as 
compared with solutions to eq. (p^ . 

Now, let us specialize to the scalar GW metric eq. (f29]) . We take into 
account only the linear perturbations in terms of the scalar GW amplitude 

The first test mass is described by the central timelike geodesic x*(t) = 0. 
The correspondent tangent vector is u"- = (1,0,0,0). The second test mass 
is situated in the unperturbed position 2;*(0) = /' having zero unperturbed 
velocity }153j. We assume that the frame of the local observer is located 
along the central geodesic. The goal is to find the trajectory of the second 
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test mass using the geodesic deviation equation eq. (P0|) . The deviation 
vector can be written hke [153] 



N\t) = P + 5l\t) (42) 

where the variation in distance is caused by the scalar GW. Using the 
frame of the local observer, one can replace all the covariant derivatives in 
eq. (flOj) by ordinary derivatives \153\ I154j . In the lowest approximation eq. 
lO]) reduces to eq. and specializes to 

^ = -'-P^m = ^uJ^V^^^ (43) 



in the field of a scalar GW eq. (l29l) . The relevant solution to eq. (1431) 
coincides exactly with the usual electric part of the motion given by equation 
eq. dl]). As we want to identify the magnetic part of the gravitational force 
arising from a scalar GW, all the terms in eq. (|40p have to be considered. 
Since is of the order of <I>, the third term of eq. (I40p is of the order 
of and can be neglected. Working out the derivatives of the curvature 
tensor and substituting them into equation eq. (|40p specialized in the field 
of a scalar GW eq. ([29]) . the accurate equations of motions read 

= \uH^^ef - \uHH\h5^^ + \k^)5l^e^:l (44) 

In this equation, which clearly exhibits two contributions, the second term 
is responsible for the magnetic component of motion and can be interpreted 
as the gravitational analogue of the magnetic part of the Lorentz force (see 
also the analogy for ordinary tensor waves in previous Subsection). 



2.4 Variation of distances between test masses and response 
of interferometers for the scalar magnetic component 

It was already recalled that the previous descriptions in the frame of the lo- 
cal observer are as close as possible to the description of laboratory physics. 
As all the questions concerning test masses positions have been analysed, 
now it is possible discussing the variation of distances. We are interested 
in the distance between the central particle, located at coordinate origin, 
and the particle located, on average, at some position (Zi, hih)- This model 
represents the situation of the beam - splitter and one mirror of an inter- 
ferometer |153l I160j . In the frame of the local observer the line element is 
given by equation |154j 
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ds^ = -{dx^f + Sijdx'dx^ + (45) 
which gives the Gahleian distance 

d{t) = x/2;2 + y2 + ^2 + 0{[^l{ujl)f). (46) 

Eq. (I46p is accurate for terms of the order of and inclusive, 
while the terms quadratic in $ are neglected. Putting 

X = li + 5x 

y = l2 + 6y (47) 
z = I3 + 6z, 

we get ll53j 

d{t) = l + j{h6x + l25y + h5z) (48) 

and, using the time dependent positions eq. ([33]) . the distance d{t) is ob- 
tained with the required approximation (i.e. ul <C 1) 

d{t) = ^ + ^(^? - ^i)^M) - ll^kill - /i)^ {u:t - I) . (49) 

Clearly, the first correction to / is due to the electric contribution, while 
the second correction to I is due to the magnetic contribution. 

Now, let us compute the response of a laser interferometer. To compute 
the response function for an arbitrary propagating direction of the scalar 
GW one recalls that the arms of the interferometer are in the it and if 
directions, while the x,y,z frame is adapted to the propagating scalar GW. 
Then, once again, the spatial rotation of the coordinate eq. (fT2]) has to be 
performed. 

In this way the scalar GW is propagating from an arbitrary direction 1^ 
to the interferometer (see Figure 2). 

At this point, one recalls that the response function is given by 

Sd{t) = du{t) - d^{t), (50) 

where du{t) and dy{t) are the distances in the u and v direction, and, 
using equations eq. (HOl) . eq. (fT2]l . and eq. (f50]l it is 



5d{t) = -^{t)l sin^ e cos2(t>+<^{t)ujl^^cos 9 



1 + sin^ ( 



+ sin^ 6 sin 2(j) (cos 
(51) 
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Figure 5: the angular dependence of the total response function of the LIGO 
interferometer to the magnetic component of a SGW for / = 8000 Hz 

In this equation the first term is due to the electric contribution, while the 
second term is due to the magnetic contribution |153j . The function 

(cos — sin 0) I (52) 

represents the so-called angular pattern [153j of interferometers for the mag- 
netic contribution. The frequency-dependence in this angular pattern ren- 
ders the magnetic component important in the high-frequency portion of 
the interferometers sensitivity band. Its value is shown in Figure 5 for the 
LIGO interferometer for the frequency / = 8000 Hz, which falls in such a 
high-frequency portion. 



ujI- cos ( 
4 



1 + sin^ 9 ^ 2 ^ 

+ sm 6 sm 2(p 
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Again, from the Figure, it looks clear that if one neglects the magnetic 
contribution, approximately 15% of currently observable signal could, in 
principle, be lost [153] . 

3 Conclusion remarks 

After extensively reviewing general relativistic gravitomagnetism, both his- 
torically and phenomenologically, the so-called magnetic components of GWs 
have been reviewed in detail. Such components have to be taken into ac- 
count in the context of the total response functions of interferometers for 
GWs propagating from arbitrary directions. Following the more recent ap- 
proaches of this important issue, the analysis of such magnetic components 
has been reviewed in both of standard GTR and Scalar Tensor Gravity. 
Thus, it has been shown in detail that such a magnetic component becomes 
particularly important in the high-frequency portion of the range of ground 
based interferometers for GWs which arises from the two different theories 
of gravity, the reviewed results have shown that if one neglects the magnetic 
contribution to the gravitational field of a GW, approximately 15% of the 
potential observable signal could, in principle, be lost. 
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